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DISCUSSION OF PROGRESS REPORT ON STUDIES 
ON THE DESIGN OF STABLE CHANNELS 
BY THE BUREAU OF RECLAMATION 
PROCEEDINGS-SEPARATE NO. 280 


SERGE LELIAVSKY, ! M. ASCE—Professor Lane has produced a 
very complete summary of work on the problem of Scour in Canals to 
which only minor additions are possible. For instance, the fact that in 
Professor White’s tests on scour, the angle of repose was used in a 
more general sense than in the paper under review, may possibly widen 
the general outlook on the problems dealt with therein. In fact, whilst 
Professor Lane introduces the angle-of-repose concept starting with 
the side-slopes, Professor White begins with the -horizontal bed-surface 
and demonstrates experimentally that the scour criterion for such a 
plane is capable of being expressed in terms of this same angle. 

On the other hand, it may also be interesting to recall, that a paving, 
consisting of a coarser grained material than the average which Profes- 
sor Lane found to have naturally formed on the bed surface on the San 
Luis Valley canals, has been obtained, earlier, experimentally, in Pro- 
fessor Rouse’s well known tests.2 The result is perhaps important as 
a point in the argument in favour of the effective size of the grain being 
larger than the median. It is, consequently, desirable to continue these 
experiments, in order to ascertain the reason of the sudden change in 
conditions, when the size of the grain is reduced to less than 5mm. 

The writer wishes to emphasize that these are very minor remarks 
only, as compared with the splendid piece of research work presented 
in Professor Lane’s paper. What makes it even more valuable is the 
fact that it constitutes the first genuine attempt to rationalize a design 
problem, which heretofore has been solved by empirical (or semi- 
theoretical) equations; and therefore, even if some among Professor 
Lane’s solutions are not final, his paper opens, nevertheless, the door 
for further widespread developments of the rational canal-design meth- 
od. 
It is this particular aspect of Professor Lane’s paper which calls for 
especial attention, since, in the writer’s opinion, the paragraph entitled 
“Further Studies”, is the only part of the paper upon which opinions 
might possibly differ, for the program sketched out therein lays possibly 
too much stress on the conventional approach to the tractive-force pro- 
blem and canal designers’ may perhaps disagree with the principle of 


1. Civ. and Hydr. Engr., Cairo, Egypt. 

2. C. M. White, “The Equilibrium of Grains on the Bed of a Stream,” 
Proc. Royal Soc., Series A, Vol. 174, Experiments 7 and 9, Table 2, 
page 331. 

3. H. Rouse, “Criteria for Similarity in the Transportation of Mate- 
rials,” Proc. Hydr. Conf., Univ. of Iowa, Bull. 20, 1940. 
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canalising the bulk of the relevant research into one single channel. 

In this connection it will be realized that the use of the tractive force 
in Professor Lane’s method is determined by the following fundamental 
points: (a) the distribution and local intensities of the tractive forces 
are calculated from either isovels, velocity equations or membrane 
analogies which, presumably, are, or are meant to represent, average 
in time values, (b) the tractive forces are assumed parallel to the axis, 
(c) the “lift”, i.e. the kinetic force normal to the surface of the channel 
is disregarded, and (d) the loci of permanent local accelerations are 
ignored. In any case if the foregoing interpretation of Professor Lane’s 
method is wrong, or not accurate, many readers will, certainly, be 
grateful for additional explanations on the subject. 

In so far as they reflect on the program of future studies, these four 
postulates will next be considered seriatim. 

In regard to (a) Professor Lane has shown in an earlier publication 
that the capacity of the water to dislodge sediment grains does not de- 
pend on the average flow conditions, nor on the maximum values of the 
relevant parameters, but on the whole range of the turbulent pulsations, 
as determined by the Gaussian frequency curve; which is to be compound- 
ed in such computations, with the frequency curve of the grain sizes. 

Since there is no reference to this point in Professor Lane’s new 
paper, it is not quite clear whether the curves in his diagrams refer to 
averages or to maximum values of the fluctuating tractive force at a 
point; or as to that, to some “effective” parametrical value thereof, 
which is chosen, within the fluctuating range, in such a manner as to 
assess correctly its true quantitative scour action. 

Had the fluctuations of the local tractive forces at the individual 
points of the perimeter, been insignificant, or had they been rigidly cor- 
related with the respective averages, the point of this argument would 
not have mattered very much, for the average might have, then, been 
chosen as the significant parameter; this, however, does not seem to 
be the case. For instance Professor White concludes that in certain 
cases the maximum tractive force is 2.25 times the average, whilst 
in others it is 4.00 times the average. Should the range of variation 
be indeed as high as that, the quantitative concept of a local tractive 
force would appear illusory, unless its relation to that range were fully 
and explicitly defined. The point needs obviously to be investigated 
more fully. 

It is, therefore, tentatively suggested to include in the program of 
future studies the problem of the pulsations of the drag force, in so far 
as these pulsations affect the stability of the grains forming the channel. 
The method of this study is, of course, left to the officers in charge of 
the investigation, but the use of the ultra-microscope, as in Fage’s ex- 
periments, ® is mentioned, for information. 

Turning, next, to point (b), it will be remembered that in spite of all 
precautions taken by Professor Vanoni to produce parallel flow ina 


4. E. W. Lane and A. A. Kalinske, “The Relation of Suspended to Bed 
Material in Rivers”, Trans. Amer. Geoph. Un., 1939, p. 637. 

5. Loc. cit., pp. 328 and 333. 

6. Fage, Rep. Aero. Res. Comm., London, 1933. 


522-2 


straight flume, sediment transportation in his experiments? neverthe- 
less, took place in the form of bands, which were symptomatic of the 
existence and the sediment-carrying capacity of well-defined, permanent 
and regularly arranged cross-currents. This case is one instance only, 
among a large variety of facts showing that the problem in hand is bas- 
ically three-dimensional, viz., that, depending on the specific circum- 
stances of the case, diagonal tractive force might (or might not) be more 
important than parallel forces. Hence, the number and shape of second- 
ary circuits for any given type section is another item of research, 
which may, perhaps, be proposed for the program of future studies; 
inclusive of the theoretical approach, as initiated by Prandt18 and of the 
experimental methods such for instance as applied by Lossievsky.9 In 
fact, in studying this problem we should not ignore the fact that Profes- 
sor White observed!0 that the granular bed in a diverging flume was 
extremely active, in spite of a negligible mean drag. Disregarding the 
non-parallelism of the flow does not seem, therefore, the attitude to be 
adopted in working out the program for the new proposed studies con- 
cerning the effect of the drag on the granular channel materials. 

It is, however, the point (c) which deserves particular attention in 
connection with these studies, for no conclusive proof has thus far been 
produced that “lift” was to be disregarded, and scour analysis based 
exclusively on tangential forces. 

In fact, “lift”, or the kinetic hydraulic force normal to the surface of 
the channel, appears, as the center of interest, in a number of papers 
on sediment movement and allied problems, and, although some among 
the modern writers on the subject11 quote Professor White as having 
demonstrated that lift is non-existent, or, possibly, negligible, the meth- 
od of this demonstration is subject to considerable doubt. For, had there 
been no vertical component transmitted by the string to the experimental 
wax grain which was attached to it and placed, in this experiment, on the 
bottom of the experimental flume, this grain must have risen to the sur- 
face of the flowing water, due to the action of the buoyancy alone, quite 
apart from the effect of lift. 

What actually occurred was than any tendency of this grain to rise 
was prevented by the impact of the water on the string—a well-known 
design fallacy, common to the numerous “tachometers”, which, in the 
seventeenth century, were used (or misused) in Italy, and impeded to 
so great an extent the advance of the hydraulic science of the period.12 


7. V. A. Vanoni, “Transportation of Suspended Sediment by Water”, 
Trans. ASCE., Vol. 111, p.67. 
8. L. Prandtl, “Uber die ausgebildete Turbulenz”, Proc. 2nd. Int. Congr. 
for Appl. Mech., Zurich, 1927, p.62. 
9. Cf. M. Welikanov, “Flussbettbewegungen und Geschiebefuhrung”, 5th 
Hydr. Confer. of Baltic States, Helsinki, 1936. 
10. Ibid. p. 333. 
11. Cf. Kalinske, “Movement of Sediment as Bed Load in Rivers”, Trans. 
Am. Geoph. Un., Vol. 28, No. 4, p. 616. 
12. Serge Leliavsky Bey, “Historic Development of the Theory of the 
Flow of Water in Canals and Rivers”, The Engineer, London, 
April 13, 1951, p.467. 
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The fact, therefore, that Professor White’s grain did not rise to the 
surface of the flowing water, cannot be taken to provide evidence that 
lift is negligible; particularly since, according to his own statement,13 
Fage’s unpublished experiments show that “lift” is at least as large as 
drag. 

Various mechanical “models” have been suggested in order to explain 
the action of the lift; for instance, according to Einstein!4 and Rubey 
lift is a matter of pulsations, whilst Jeffreys!6 conceives it as a com- 
ponent of temporal average forces. The trouble about all these “models” 
is that whilst obvious to the scientist, they do not appear convincing to 
a practising designer. A force parallel to the direction of the flow is 
evident, but, unless its nature is clearly demonstrated, a force normal 
to the flow, is an abstraction. 

From this standpoint an elementary but obvious “model” explaining 
the kinetic “lift” is desirable. 

Taking, therefore, as datum the horizontal (or nearly horizontal) 
plane AB (see Fig. 1 on the annexed Plate 1), representing the average 
surface of a granular channel bed, the relative position of an individual 
grain belonging to, or forming part of, this surface, may be defined as 
follows: highest, when the particle is above the plane (as in a), and 
lowest, when it is entirely below the plane (as inb). In the average 
position of the particle, the plane must, therefore, strike its center (as 
in c). 

This average position can, consequently, be taken as a starting point 
for a preliminary “statistical” examination of the problem. The advan- 
tage of this approach is that we can, thus, draw a parallel between the 
phenomena of scour and sphere-drag, and obtain a general idea about 
the former, from the abundant information available about the latter; 
for instance, in Fig. 2 (see same Plate), which is based on this analogy, 
we discriminate between three possible flow patterns near and around a 
sediment grain, which show: (a) that the drag and the lift, cannot be re- 
presented by one and the same parameter applicable to all the phases of 
flow and to all the shapes and arrangements of grain; (b) that changes 
from one flow pattern to another are likely to occur suddenly (as they 
sometimes do in experiments with sphere-drag), which furnishes a 
plausible explanation of the too often disregarded fact that the intensity 
of the scouring action does not always develop gradually, as must have 
been the case had the traditional drag formulae been rationally correct. 
But, as observed in numerous experiments, scouring action undergoes 
abrupt and apparently anomalous alterations, which however can only 
be due to sudden changes in the flow patterns for the bulk (in the statis- 
tical sense of the term) of the grains, (c) that the vertical component of 
the kinetic force applied to the particle is directed upwards, correspond- 
ing to the vacuum pressures observed in the sphere-drag experiments. 


13. Loc. cit. p. 333. 

14. H. A. Einstein, “Formulas for the Transportation of Bed Load”, 
Proc. ASCE, March 1941, p.351. 

15. W. W. Rubey, “The Force Required to Move Particles on a Stream 
Bed”, Prof. Mem., U. S. Geol. Survey, No. 189-E, 1938. 

16. H. Jeffrey, Proc. Camb. Phil. Soc., 1929, 25, 272. 
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This vertical component we describe as “lift”. 

Better to understand the argument of this last point, attention is 
directed to Fig. 3 in the same Plate, which represents an experiment 
published by the writer as early as 1914.17 The model of an air-ship 
hanger was tested in a wind tunnel, in order to find the wind forces for 
its structural design. The point of the experiment was, that, in spite 
of the wind direction being inclined to the horizontal (10° subtangent 
angle, as shown in the diagram), the resultant wind force applied to the 
building was directed upwards, and not downwards, as the standard 
specifications would have it. 

This result is in full agreement with the well known fact that, during 
severe storms and hurricanes, roofs of buildings are never pressed 
downwards, but are frequently lifted from their original positions and 
transported to surprisingly long distances therefrom. 

Should another obvious and elementary application of the same prin- 
ciple be deemed necessary, or desirable, as another illustration of the 
lift model, we may cite the common vaporizer, shown in Fig. 4. It 
consists of a combination of two tubes, one being vertical and the other 
horizontal, which are held together, by some means or other, close to 
the point where their axes meet. By blowing through the upper pipe, a 
horizontal air current is created, which suffices to cause the liquid to 
rise vertically to the top of the lower pipe, and be, then, diffused in the 
form of a cloud. Here we have, therefore, another example of a hori- 
zontal flow, causing a vertical movement of a heavier substance. 

These simple examples are believed to supply an irrefutable proof of 
the very real existence of the lift force; and moreover, since the hori- 
zontal component of the wind pressure acting on the shed in Fig. 3 is 
obviously a function of the height of the vertical wall, whilst the upward 
component of the force applied to the roof depends evidently on the 
shape of the latter, it is clear that these two forces cannot be represent- 
ed by one and the same parameter. It follows that the lift problem must 
be envisaged as a self-contained item of research, correlated but not 
identified with the drag problem. 

In view of the foregoing consideration it may possibly be recommend- 
ed to consider the lift problem as a desirable item of the program for 
future research, in addition to, and parallel with, the drag problem. 

With reference to the last point (see point d in page 2) attention is 
called to the fact that sediment transportation is a composite process 
consisting as it does of three distinct phases. These phases are: (I) 
The initial period when the particle is picked up and moves generally 
at a lesser velocity than the water; (II) the transportation period proper, 
when both velocities are of the same order of magnitude, and (III) the 
last stage during which the particle is dropped. 

That the tractive force suffices to explain, and define, sediment 
behaviour in the first period, is a debatable but possibly acceptable 
proposition. But that it should also be taken as a significant parameter 
for the second phase, is more difficult to agree with. Since (with the 
possible exception of du-Boys’ concept of a solid mass of granular bed 


17. Magazine of the Institute of Engineers of Ways of Communication of 
Emperor Alexander I, St. Petersburg, 1914 (in Russian). 
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material shifting bodily downstream below the physical bed of the chan- 
nel) the first phase cannot be conceived without at least some reference 
to the second, it seems that the specific factors conditioning this second 
stage are also worthy of attention. 

Hence, in addition to the erosive (in the narrow sense of the term) 
power of the stream, its “capacity to carry sediment of a certain size 
range” is, possibly, still another indispensable section of a consistent 
program of silt-stable-canal research. 

In the writer’s opinion, (thus-far unpublished), apart from the canal 
section as a whole, further research should concentrate on the local 
sediment carrying capacities of individual points of this section, in 
correlation with its typical pattern of secondary circulation; for the 
concept of uniform flow as presented in the elementary textbook, has 
outlived its time; and even if we deal with temporal averages of a gen- 
erally pulsating current, secondary circulation must of necessity cause 
local foci of permanent accelerations and decelerations, even though the 
average velocity, taken for the whole section, were constant over the 
length of the investigated reach of the canal. 

From considerations bearing on the continuity equation, it will then 
be found that acceleration capable of a greater sediment carrying capaci- 
ty must be correlated with convergent currents, whilst deceleration, and 
therefore sediment deposits, must coincide with divergent currents. 
Hence, stable canal form is very probably a function of the cross-current 
pattern, and since Lossievsky! found that this pattern changes with the 
bed-to-width ratio, we have at least a glimpse of an idea as to the me- 
chanical processes which lie behind the tendency of a stream to exca- 
vate a channel of a certain, given geometrical form. 

This explanation, preliminary as it may be, may nevertheless be 
considered as an alternative to the concept of “incipient motion” profile, 
described in the paper under review, and deserves to be investigated in 
parallel with it. 

A. NIZERY?9 and G. BRAUDEAU29—pProfessor Lane’s study makes 
a clear distinction between the two basic problems: (1) What is the 
necessary hydrodynamic effort to put in motion a grain of a given mate- 
rial? (2) What form of section must be chosen in order that the forces 
applied on the material be distributed as uniformly as possibie? 


I - The first problem may be studied in a very wide canal with horizon- 
tal bottom. Thus the gravity will have no tangential component and 
the tractive force & , the value of which is well known in this case, 
will be the only one acting on the grain, parallel to the wall. The 
grain equilibrium will be realized if the resultant of the tangential 
force and of the gravity is situated inside the friction cone defined 
by a certain angle ¥ . The tangential force is proportional to T 
and to the grain surface, and consequently to d2. The weight is pro- 
portional tod3 and YY, (density of the material under water). The 
equilibrium condition is therefore: 


18. Loc. cit. 

19. Directeur Adjoint des Etudes et Recherches d’ELECTRICITE DE 
FRANCE. 

20. Ingénieur au LABORATOIRE NATIONAL D’HYDRAULIQUE - Chatou - 
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When the grain sizes are not too different, one can consider for 
all the grains only one angle ¥ , equal to the natural bank angle of 
the materials (under water). The law found in this case is again 
similar in form to that deduced from the MEYER PETER formula, 
which is currently used in Europe to study the problems relating to 
bed load transportation. It shows that the critical tractive force 
T. is proportional to the diameter and relative density of the 
grains. 

We often have had the opportunity of applying this semi-empiri- 
cal law, which results from a very large amount of experimental 
work. We always found a very good agreement with experimental 
results, when the conditions of application for which this law had 
been established were suitably satisfied. When however the grain 
sizes are too different, or if cohesion comes into play, there is no 
more agreement and it then becomes necessary to resort in each 
case to direct measurement of U,¢. 


As shown by Professor Lane, the best form of section will be that 
for which the tractive force at each point is equal to the critical 
force for the materials with which the walls are built up. It is 
such a section which will allow the greatest mean speed for a given 
rate of flow, and will be therefore the most economical as regards 
the excavation work. 

If the tractive force distribution were known for a determined 
form of section, it would be possible to express that the force equi- 
librium is satisfied at each point, in a form similar to that indicated 
above. It would then be possible to derive the equation of the section 
which complies with this condition, and also in a general manner 
with the five conditions indicated by Professor Lane. 

But the tractive force distribution depends on the form of section 
and it follows that the equation found would be theoretically an in- 
tegral equation. 

One is therefore compelled to make hypotheses and to check if 
these are confirmed experimentally. 

These questions are now being analysed at the Laboratoire 
National D’Hydraulique, at Chatou, on a canal with a length of 45 m 
and a width of 3,50m. The form and size of the section, as well as 
the longitudinal slope of this canal can be easily varied. As we are 
only considering up to now the forms of section, we use as a mate- 
rial an uniform sand, the grains of which have very nearly 3 mm of 
diameter. 

For a given section form and size, we increase the slope, while 
modifying the rate of flow at the same time, so as to obtain the 
uniform regime, till a bed load transportation clearly begins to 
occur at any point of the section. 

For each regime, velocity measurements are effected at many 
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points. We have thus been able to obtain detailed velocity distribu- 
tions for different types of sections, especially for trapezoidal ones. 

From these velocity distributions the tractive force at each point 
has been calculated as indicated by J. B. Leighly and Professor 
Lane. 

We have found that the experimental results did not agree with 
those obtained in the above manner: the bed load transportation did 
not begin where the calculation indicated a maximum tangential 
force. Perhaps this phenomenon is not alien to the existence of 
secondary flows above the section corners, as shown by the strong 
instability of the velocity measured in these zones. 

On the other hand, we have calculated a section in which the 
tangential forces would be uniformly distributed, on the hypothesis 
that the tractive force at any point would be simply: 


(1) 


with Z = depth of the considered point below the surface, and 
S = slope of the water line. 
It is easy to find for such a section parametrical equations of 
the form: 


x = a cos 


sing / 1- sin’ y aa (2) 


(3) 


kyé 


Here P = weight of a grain 
K = coefficient of proportionality of the force applied by the 
water on a grain to the tractive force at the same point. 
a = maximum depth. This is correlated with the slope s by 
equation (3). 
It is further necessary to add a friction equation correlating a 
and S with the rate of flow.21 
We have made experiments with such a section, in the same 
manner as for those considered above, and found that the bed load 
transportation begins almost simultaneously at all points. This 
would indicate therefore that the force distribution is fairly well 
represented by eq. (1), at least in a section, the height/width ratio 
of which is below a certain value, which, as a matter of fact, we 
have not determined. 


III - It seems that one can expect from such studies more than the mere 
development of a rational method of calculation for canals with al- 
luvial bottom. 


21. Another hypothesis, made by KOECHLIN about the manner in which 
the water acts on the grains, leads to a section, the equation of which 
has the form Z =a . sin (tg¥ z). It is found that the sections cal- 
culated by both methods for the same regime are rather similar, 
though of different analytical expressions. 
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The application to natural watercourses of the results so obtain- 
ed may be fruitful. If we consider a river meandering in a high- 
water bed built up by its own alluvia, observation suggests that the 
hydrological regimen assign to the low-water bed a certain slope. 
We have thus four variables, which do not appear to be independent, 
since the rate of flow values, for example, must determine the val- 
ues of the three others. The mere knowledge of the friction law and 
of a specific bed load transportation law22 does not enable one to 
account for this occurence. The width must appear in the lacking 
relation, and it seems that this cannot be other than that which ex- 
presses a certain (uniform, or non-uniform since curvature effects 
arise) distribution of the tractive forces, necessary for the trans- 
verse profile equilibrium. 

The problem is of the same nature, although naturally more 
complex, owing to the fact that on one hand the rates of flow vary, 
and that on the other hand we consider an equilibrium with bed load 
transportation. 

We believe that considerations similar to those used to study the 
equilibrium without bed load transportation under permanent regime 
may help to solve this important problem. 

JOSE S. GANDOLFO23—The author’s contribution describes new and 
improved methods in the design of unlined canals by means of the trac- 
tive force concept. 

It gives the results of the intensive experimental works conducted at 
the Bureau of Reclamation under the author’s supervision. 

It also mentions observations, investigations and technical experiences 
of other countries, especially India. 

The author’s methods for designing stable canals in coarse and fine 
noncohesive material represent a practical arid valuable advance in 
hydraulic projects. 

There are many aspects of these problems which need further studies 
and investigations. Among them, the writer considers important to com- 
ment some ideas presented under the author’s heading: “Canals Carry- 
ing Heavy Sediment Loads”. The following discussion is concerned pri- 
marily with this part of the author’s paper. 

These canals are often in connection with natural streams and un- 
steady flow will be expected. If unlined, the bed is of the same nature 
as the heavy load material transported. 

In this special case we cannot speak of stable channels because each 
flow has a definite mean rate of bed load transport and slope. So, in 
periods of lower flow, ordinary flood deposition and decreasing of gen- 
eral slope will be expected. In periods of great flood there are erosion 
and increasing of slope and rugosity. Normal flow can not remove all 
sizes of coarse bed load and only transports the small particles not 
protected in spaces between coarser particles. 


22. We call “specific bed load transportation law” the law relating the 
solid flow per meter of width to the hydraulic conditions relative to 
an infinite width. 

23. Head of Hydraulic Department and Professor of the Universidad 
Nacional de Eva Perén (Ex-La Plata) (Argentina). 
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Since the canal configuration is changing continuously, it is neces- 
sary to choose a representative flow for correct design. 

Hydraulic needs will be satisfied by maximum flood flow canal capaci- 
ty, but alluvial requirements are connected with the necessity of estab- 
lishing the prevailing flow of coarse transported material. 

The product of the flow frequency and the respective canal coarse 
sediment transported give the frequency curve of the coarse sediment 
transported; this shows a maximum ¢G' for a determined flow Q'. 

The narrow band of flow above and below Q' carry the greatest 
amount of heavy material load. Then, the canal slope and end bed 
Contigueatten will agree well with this value of Q' (in Fig. A, Q' 

650 m°/s). F. Schaffernak,24 named it “bed generative flow”. 

In the dike improvement design of the San Juan River (Argentine) the 

writer found that the mean annual heavy material transport was: 


G' m = 103 kg/s; corresponds to a flow: 
Q' m = 538 m/s, but the mean annual flow is: 


= 521 m3/s 


The bed generative flow concept gave the generative coarse sediment 
transport as G' = 175 kg/s, and the corresponding flow as 


Q' = 650 m3/s 2 


The flood flows between 450 and 850 m°/s transports 75% of the mean 
annual load; also, G' = 1,7 G' m- 

The Rhine River at Brugg Bridge, above Constance Lake26 showed a 
mean annual transport during 1936 of 


G' m = 7,3 ke/s corresponding to a flow of 
Q'm = 503 m/s ; 


but the mean annual flow during 1936 was 
Qm = 457 m3/s 


The bed generative flow concept gave the generative coarse a, 

load as G' = 16,1 kg/s and the coprenponans flow as Q' = 659 m°/s. 
The flood flows between 520 and 940 m3/s transported 62% of the 

1936 coarse sediment load; also G' was equal to 2,2 G' m. 


24. F. Schaffernak: Die Theorie des Geschiebebetriebes und ihre Anwen- 
dung. Zeitschrift des Osterreichischen Ingenieur and Architeken- 
Vereines, Wien 1916, Nr.68- Neue Grundlagen fiir die Berechnung 
der Geschiebe fiihrung in Flusslaufen. Verlag:Franz Deutike.Leipzig. 
Wien-1922- 

25. José S. Gandolfo-Estimacién de la evoluci6n del alveo de los rios. 
Publicaci6n de la Fac.de Ciencias Fisicomatematicas, Serie Segunda, 
3, NO 128, pag 245 y siguientes, Separata pag. 33 y siguientes, Uni- 
versidad Nacional de La Plata (actualmente Eva Per6én)-Argentina, 
1940. 

26. F. Nesper-Die Internationale Rheinregulierung von der Illmiindung 
bis zum Bodensee-Schweizerische Bauzeitung, Zurich 1937- Band 


110, Nr.13. 
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27. E. Meyer-Peter, H. Favre u R. Miiller.Beitrag zur Berechnung der 


The design for International Rhine regularization from Illmiindung to 
Bodensee was computed with G’ ,,; in the writer’s opinion it is an in- 
sufficient dynamic design basis for rivers.27 

If the flood-flow frequency would be constant the mean annual load 
basis will be correct; on the contrary, heavy material bed transport 

increases as flood flows decrease in frequency; the values above are 
illustrative. 

The differences between G'm and G' will be greater with the higher 
irregularity of the number and distribution of flood periods; also with 
the more variable intensities and duration of the flood waves. 

Then, these differences will be greater in the improvement works 
of torrential rivers. 

The mean annual bed load, if the other conditions are the same, shows 
a longitudinal profile of less slope and a characteristic cross section 
not deep enough. 

Otherwise, in the mean annual coarse sediment bed load is included 
all the transported material, corresponding to lower and greater flood 
flow. For that reason, this value is highly variable. 

The generative bed load flow is a value little influenced by flood flows 
weak in solid transported material, even if the frequency is increased. 

This value is practically not modified by the extraordinary floods of 
low frequency, though they transport a great lot of solid bed material. 

The generative bed-load flow is confined within the narrow flood- 
flow zone, which because of its frequency and capacity for transporting 
solid bed material moves the greater volume of coarse sediment. 

Generative bed load flow is much less variable than mean annual 
flow. 

In other consecutive canal stretches generative bed load flow can 
change through modification of the granulometric curve of the bed mate- 
rial or slope. This is also true for a varying flow when a tributary 
canal reaches the beginning of the following stretch or where there is 
a diversion structure. 

Problems of designing canals with transport conditions that vary in 
consecutive stretches can be attempted as follows: 

Let: generative bed load flow = Q', coarse material bed transport = G', 

mean diameter = d, specific gravity = ? , dimensions of the geo- 

metrical and hydraulic canals be B (width), ‘? (h) (height), and i 

(slope); and let k = rugosity. 

Use subscript 1 in the first stretch of the canal and the heavy bed 
material transport function: 


If G'1 is known and if ij is the actual slope or another slope conven- 
ient for design, then solve the function by steps varying ~ (hy). Imme- 
diately compute: By, = (Q'1, f (hi) , ky ): 

From the variation of both iy and f (hy) in the G' 1 equation, the best 
B, value can be chosen. 


Geschiebefiirung und der Normalprofilbreite von Gebirgsfliissen, 
Schweizerische Bauzeitung, Zurich, 1935, Band 105,Nr.10- 
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In the second stretch of the canal, with subscript 2, the continuity 
law of transport of heavy material for the width of the canal is: 


ay - 
\P 
in which | represents the part of heavy transported bed 


material and the rest being crumbled and incorporated as suspension 
flow; when alluvial material is added at the beginning of this stretch G2 
has a positive value; but when there is a diversion structure or there is 
extraction of the coarse material, or the generative bed flow diminishes, 
Gg has a negative value. 

We must investigate, in all these cases, the value of Q' 9, but, if 
Gp = 0, it is accurate enough toestimate Q'yg = ~~ Q', . 

If there is a discontinuity in the heavy solid material transported in 
any stretch of the canal (G, 40), the corresponding G',; and dy values, 
will be the basis for solving G' in the consecutive stretches of the canal. 

In the second canal stretch, the slope ig can be solved: 

' 
4 2 


The second member is known by giving values to 4 (ho) = 
However: 


(i) = ft, Fai | 5 


by giving values to i, f2(i) can be solved and the curve for|i,f. (i) 
drawn. 

In this diagram the ordinate f = (45) gives the slope ig. 
Finally: 3, = p “kp 
The best By value can be chosen from the variation of f (ho) ; solve 
again forfg (i) and drawn the new diagram [i, fo (i)] 
Immediately other values of ig and Bg can be computed. 

After the dimensions of the cross section of the canal for heavy mate- 
rial transport have been computed the hydraulic capacity for maximum 
project flow must be computed. 

The writer, in the contribution25 gives examples of several charac- 
teristic cross section applying Meyer-Peter, Favre, Miiller’s formula 
for heavy solid transport. The procedure is general, and can be applied 
with any other heavy solid transport formula for computing G'. The 
Einstein formula2®8 js the most advanced, but of laborious application. 
In the writer’s project of an improved dike for the San Juan River, 


28. H. A. Einstein—The bed load function for sediment transportation 
in open channel flows, U.S. Department of Agriculture, Technical 
Bulletin No. 1026, 1950. 
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Argentine, 29 the heavy solid material transport continuity, in the con- 
secutive stretches was determined using the Meyer-Peter, Favre, 
Miuller’s formula, the most advanced at that time. 
. The successive stretches have conjugated slopes; the lineal function 
is: 
ig =a iy - b 
in which a and b depend on values of Q'1, Q' 9, dy, dy,F 4, 2, By, Be- 

The writer considers author’s contribution a great advance in the 
design of stable channels and hopes that soon other new contributions 
can be added to the great amount of themes that are mentioned by the 
author under the heading “Future Studies”. 

N. K. BOSE20—The idea of introducing tractive force analysis in the 
design of stable channel in alluvium is not a new one, but the lines of 
attack followed by the authors are certainly interesting. The basic idea 
“that in all canals having the same ration of B to D and the same side 
slopes, the tractive force distribution will be similar” requires more 
elaboration and proving. The authors have tried by various methods to 
determine the tractive force distribution along the wetted perimeter of 
a canal section from various standpoints. The use of the membrane 
analogy method appears to be capable of giving satisfactory solutions 
and may help in the design at a later stage. 

However, the main difficulty experienced so far by the different re- 
search workers in developing a theory of regime channels has been 
mainly in correlating silt or sediment movement with the other hydro- 
dynamic characteristics of the channel flow. Whether we use tractive 
force or any other hydraulic constant of a canal, such as, mean velocity 
“Vm”, hydraulic mean radius “R” wetted perimeter “P”, area of the 
waterway “A” or depth of water “D”, the greatest difficulties has been 
to bring in a suitable characteristic of the sediment that the canal car- 
ries on the bed or in suspension. The authors’ analyses also donot show 
any advance on this aspect of the question. Their assumptions of taking 
K25 or Kso as representative diameters of the sediment have not been 
justified. 

The extensive work that was carried out by the research workers in 
the Irrigation Research Institute, Punjab, India, during the years 1932- 
42 have shown conclusively that the sediment that a channel carries can 
have great influence in determining the ultimate section and slope to 
which a channel will settle down. This fact was most probably in the 
mind of Kutter and Ganguillet when they tried to modify the Chezy for- 
mula by introducing in it a variable factor “n” which depends on the so- 
called rugosity of the channel. About 50 years ago Mr. Kennedy of the 
Punjab Irrigation, introduced the idea of silt and turbulence explicitly 


29. José S. Gandolfo—Estudio de la evoluci6én fluvial que determina el 
endicamiento del Rio San Juan-Publicaciones Especiales NO 2- 
Facultad de Ciencias F{sicomatem&ticas de la Universidad Nacional 
de La Plata (actualmente Eva Perén) Argentina- Enero de 1940, 

NO 126, pag. 97. 

30. Director, River Research Institute, West Bengal, on the progress 
report of studies on the design of Stable Channels by the Bureau of 
Reclamation. 
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into the subject and deduced his well known relationship - 
Vo = KD" 


More recently Mr. Lacey of the United Provinces, India, had again re- 
vived the subject and proposed in his two papers before the Institute of 
Civil Engineers, London, a number of relationship involving some of 
the usual hydraulic constants with a quantity known as silt factor “f”. 

It is well known that the main difficulty that the previous investigators 
had in dealing with this problem of regime channel was in defining silt. 
Mr. Kennedy in his siltometer proposed an instrument, which, but for 
some minor defects would have enabled him to classify silt successfully— 
some of these defects have since been removed by Dr. Puri and the 
modified siltometer is used extensively in India. The siltometer that has 
solved the problem of classifying silt and made it possible to go ahead 
with silt investigation had been perfected by Mr. Vaidyanathan. A paper 
No. 167 on this siltometer was read by its author before the Punjab En- 
gineering Congress in 1933. 

In a paper (No. 192) before the Punjab Engineering Congress in 1936 
on silt movement and design of channels the writer had examined the 
silt data obtained from a number of regime channels in the Punjab col- 
lected under stable conditions. The distribution curves of the silt sam- 
ples were analysed and sub-divided into different types. These types 
could not however be put in absolutely cut and dried compartments. 

Very often they merge one into another. Just as a spectrum of a source 
of light runs from the ultra violet through the visible part to the infra 
red end, silt curves of a canal also run through the above types if they 
are analysed throughout the entire year. Each of these types indicates 
a certain condition of equilibrium of a canal, so the annual frequency 
distribution of these types of curves for any canal is a very reliable in- 
dication of the stage of equilibrium of that canal. Now the first problem 
was to determine what was the significance of these types, to what stage 
of the development of a canal they belong — what indication they give 
us of the stability of a canal. 

After examining these silt curves and the hydraulic data for a number 
of stable channels the writer, in another note on the design of channels 
in alluvium, Paper No. 252, the Punjab Engineering Congress, 1942, had 
derived the following relationship: 


A= 1.15 Q°85 
D= 0.39 Q:29/89-57 ... (2) 
Sx10° = 2.09 m:85/q:2!_ (3) 


where A is the area of waterway in sq. ft., Q the discharge in cu. ft. per 
second, D is the average depth in feet, S is the slope of the water surface, 
m the mean diameter in mm of the bed sediment. 

These three equations will give the design engineer all the informa- 
tion that he may possibly require to design his canals. ‘ 

In this connection it must be stated here that in all these analyses 
carried out by the writer it had not been possible to take into account 
the effect of the quantity of bed movement on the regime equations. 
Though, efforts have recently been made by workers like Einstein, Shields 
and others to introduce quantity into the regime equations, no practical 
solution has yet been obtained. 
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E. W. LANE,?1 M. ASCE—The author has been very gratified by the 
valuable information and constructive suggestions furnished by the dis- 
cussors of his paper. They have added materially to its value. 

The information supplied by Mr. Leliavsky, MASCE, regarding Pro- 
fessor White’s test is valuable, as these tests are the only hydraulic 
experiments known to the author which test the hypothesis of the rela- 
tion of the angle of repose to the transportation of coarse sediments. 

Mr. Leliavsky has listed four areas in which uncertainty exists re- 
garding the reliability of assumptions used by the author in his canal 
design procedure and makes suggestions for further studies to clear up 
these uncertainties. As Mr. Leliavsky stated, the values of shear used 
by the author were the average values with respect to time, and not the 
maximums. In our present state of knowledge the latter is not known in 
enough detail to be used at present in a design method. The author is 
very much in favor or further study to ciear up the question raised of 
the relation of the averages of the shear value to the maximum. 

Mr. Leliavsky has also raised the question of secondary currents, 
and suggested studies along this line. The possibility of secondary cur- 
rents being important is also suggested by Messrs. Nizery and Braudeau. 
This subject should also be added to the list of further studies mentioned 
by the author, since at present our knowledge of them is insufficient to 
enable them to be included in a design method. 

That lift has its part in sediment transportation is becoming more 
clearly evident. Here also our present knowledge is insufficient to in- 
clude it in a design method, and further information is therefore desira- 
ble, as Mr. Leliavsky suggests. The action producing lift does not al- 
ways cause an upward force. In coarse, non-cohesive material, due to 
the tendency of flat particles to place themselves in a shingle like posi- 
tion, the force due to the curvature of the currents, may be downward 
instead of upward.32 

As Mr. Leliavsky suggests the presence of loci of permanent local 
accelerations no doubt exist in some cases,”” and may be quite impor- 
tant generally. The author welcomes study of this matter, and any other 
aspect of sediment transportation, for the future progress in the im- 
provement of methods for the design of stable channels will be dependent 
or. the rate at which our knowledge of the laws of sediment transporta- 
tion increases. 

The discussion of the late A. Nizery and G. Braudeau regarding the 
shape of the channel of minimum excavation and width is very valuable. 
This supports the use of the principles involved in this study as the 
basis for design of stable channels for conditions where this method 
applies. The author has been making further studies along these lines. 
It is greatly to be hoped that a complete report on the studies of Messrs. 


31. Prof., Dept. of Civ. Eng., Colorado Agri. and Mech. College, Fort 
Collins, Colo. 

32. Some Observations on the Effect of Particle Shape on the Movement 
of Coarse Sediments. E. W. Lane and E. J. Carlson. Trans. Ameri- 
can Geophysical Union, Vol. 35, No. 3, June 1954. pp. 458. 

33. A New Method of Sediment Transportation. E. W. Lane. American 
Geophysical Union, Vol. 25, Pt. 6, p. 900, May 1954. 
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Nizery and Braudeau will soon be forthcoming. 

The lack of agreement reported by Messra. Nizery and Braudeau of 
the tractive force distribution computed by the Leighly method with the 
observed results is interesting. This may account for some of the scat- 
ter of the results observed when this method was investigated under the 
author’s direction, and supports the author’s decision to use the results 
of the mathematical and membrane analogy methods in his proposed de- 
sign method. The author agrees with Messrs. Nizery and Braudeau’s 
idea that secondary currents are probably at least a partial explanation. 
The suggestions of Messrs. Nizery and Braudeau that the principles of 
stable channel design may be used to advantage in the study of regime 
conditions for natural streams is believed to be fully justified. For 
some time the author has been carrying on a study of the form of natu- 
ral streams, and finds the principles derived for stable channels of great 
help in working out an analysis of this problem. 

Professor Gandolfo’s method for design, as illustrated by his discus- 
sion, shows that years ago he applied the principle which the author con- 
siders fundamental in the design of stable channels designed to carry a 
heavy sediment load, namely, that for such conditions it is necessary to 
determine the magnitude of the sediment load and design the canal to 
carry it. So far as the author has seen, this is the first application of 
this important principle to the design of an irrigation canal. 

Under the author’s definition of stable channel, it is not necessary 
that such a channel be always free from scour or deposit, as Professor 
Gandolfo’s discussion suggests. A stable channel, under the definition, 
may alternately deposit and scour, provided the scour carries away the 
deposit before it has accumulated in objectionable quantities, and the 
deposit refills the scoured areas before they scour sufficiently to be 
objectionable. 

Dr. N. K. Bose raises the question of the reliability of the assump- 
tion of the similarity of tractive force distribution in channels of simi- 
lar cross section. The justification of this assumption is that in all but 
very small channels the Reynolds number is so great that the effect of 
the viscus forces is negligable. Under these conditions the velocity 
distribution is similar in geometrically similar sections. The geomet- 
rical similarity of velocity distribution from model size to prototype 
size, where the Reynolds number is sufficiently large, is commonly as- 
sumed in making prototype predictions by means of models. If it holds 
from model to prototype size, it should hold over a wide range of canal 
sizes. 

Dr. Bose’s remarks deal principally with the case where considera- 
ble quantities of sediment are transported by the canals. The studies 
presented by the author had not progressed far enough in this area to 
go much beyond a statement of the fundamental principles that the author 
believes apply to the design of stable channels under these conditions. 
These are that the canal must be so designed that it will transport sub- 
stantially all the sediment brought into it, without producing sufficient 
tractive force on the banks and bed to cause scour. To determine the 
conditions necessary to transport the required quantity of sediment, 
use can be made of the formulas developed for the amount of sediment 
which will be transported under certain conditions. The greatest need 
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in this field is the measurement of the size and quantity of sediment 
transported in a considerable number of canals, so that the reliability 
of the various available formulae can be appraised, and if necessary 
new formulae developed. 

The author has not been able to secure a copy of the paper describ- 
ing the derivation of the formulae for canal design which Dr. Bose has 
proposed, and therefore he is unable to appraise it. However, it was 
based on extensive experience and probably a great deal of data. It is 
not unlikely that for the conditions in the canals which were used in 
developing it, it is under present conditions a more reliable method of 
design than one based on a computed quantity of sediment transporta- 
tion. As the laws of sediment transportation become more fully known 
however, the author expects the analysis based cn the computed quantity 
of sediment transportation to eventually become more reliable than Dr. 
Bose’s equations. For conditions of sediment load differing widely from 
those in the canals used as data in Dr. Bose’s study, the author believes 
that his equations and any other design method which neglects the effect 
of sediment concentration, will not always give reliable results. 

The writer is indebted to Mr. Leliavsky for informing him, in a pri- 
vate communication, that the solution mentioned in the author’s paper 
for the shape of a channel in coarse, non-cohesive material, which would 
require a minimum excavation and width, was previously worked out by 
P. Forchheimer*4 as early as 1930 and perhaps earlier. Forchheimer 
included in his analysis not only the effect of the slope of the sides but 
also the influence of the longitudinal slope of the canal. It is doubtful 
however if this latter effect would ever be sufficient to significant. 

The author has also discovered that the effect of side slopes on the 
movability of coarse, non-cohesive material involving the same basic 
idea as that presented by the author, but using velocities instead of 
tractive forces was published by A. Gostev?® in 1928. 

Since writing this paper, the author has been impressed by the im- 
portance as a stablizing factor of the deposits of fine material on the 
sides of many canals. Because of the cohesive nature of these deposits, 
the tractive forces which the sides can stand may be considerably great- 
er when they are present. This may occur naturally in the aging pro- 
cess which takes place when the flows discharged down the canal in- 
crease, as more and more land is irrigated by the canal but, may also 
be furthered by scientific manipulation of the natural silts. It may be 
possible to accomplish this even more effectively by the applicating of 
bentonite. A thorough study of these possibilities is being started at the 
Colorado A & M College, Ft. Collins, Colorado, in conjunction with an 
extensive investigation of the possibilities of the prevention of seepage 
out of canals by staunching the bottom and sides of the canal with a 
bentonite slurry. For canals in fine, non-cohesive material, particularly 


34. Hydraulik-Phillip Forchheimer—3rd Edition, 1930, p. 551, Copy- 
right by B. G. Teubver in Leipzig, 1930. Translation No. 47-9. Re- 
search Center, Waterways Experiment Station, Vicksburg, Mississippi. 
35. On Scouring Velocities A. Gostev. Irrigation News Bulletin (Russian) 
No. 1, Jan. 1928, pp. 65-69. (Translated by U. S. Bureau of Recla- 
mation.) 
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where the presence of a load of sand in the incoming water necessitates 
a high tractive force to prevent it from depositing, this method, if it can 
be successfully developed, will greatly improve the form of canal which 

can be used. 

In closing it is desirable to briefly summarize the views of the author 
regarding the subject covered by this paper as originally presented and 
as modified by his study of the subject since presentation and by the 
valuable discussions of the paper which were presented. The author be- 
lieves that the fundamental principles of stable channels, as given in 
his paper, are very simple and can be readily understood by any one 
with an elementary knowledge of hydraulics. Defining a stable channel 
as one in which neither scour nor deposit in objectionable quantities 
occurs,”” it follows that to obtain such a channel all that is necessary 
is to insure that at all points on the bed and banks of the canal the action 
of the water is insufficiently severe to cause objectionable scour but is 
sufficiently severe to prevent sediment carried by the water to depost 
in objectionable quantities. This statement gives the fundamental basis 
for the method of design proposed by the author. Corallary statements 
are (1) if the water entering the canal carries material quantities of 
sediment, in order that objectionable deposits do not occur, the inflow- 
ing sediment must be transported through the canal within a material 
loss of among and (2) that on the sides of a canal the force of gravity 
tends to move sediment particles down the sides toward the bottom, and 
this action assists the forces due to the motion of the flowing water in 
causing scour on the canal sides. He believes that the reliability of 
these statements are so evident that nearly everyone interested in this 
subject, after giving them athorough study, will agree to them. 

Since canal stability involves so largely the interaction of the water 
with the bed and banks of the canal, it seemed that a logical method of 
approach in devising a method of design for such channels would be a 
study in the nature of this interaction. It will be noticed that there are 
two phases of this interaction, a scour phase and a transportation phase. 
Since current knowledge of the science of hydraulics and fluid mechanics 
is not sufficient to completely analyze these phases, it is not now pos- 
sible to devise a perfect method of design, but by using the best knowl- 
edge at present available, a method can be devised on these principles. 
While he does not expect changes in the fundamentals of his method of 
analysis, he hopes for and expects that some of the assumptions which 
it was necessary to make because of the state of knowledge of the sci- 
ence when the methods were adopted, will be improved upon when more 
adequate knowledge is available. However, it is believed that the as- 
sumptions made in the author’s method are the best that is possible to 
make with our present knowledge. 

The design of stable channels for irrigation projects will not wait for 
the development of a perfect method of analysis. For practical purposes 
all that is required is that the method be complete enough to enable de- 
signs to be made by its use, sufficiently simple to be readily understood, 


36. The sides of the channel also must not slough down, but since a rea- 
sonable satisfactory analysis of this phase of stable channels is 
available, it has not been discussed in this paper. 
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logical enough to commend itself to the designing engineer, not require 
so much work as to make its application prohibitive, and be more relia- 
ble than other available methods. It is believed that for the conditions 
of small sediment loads, in its present state of development, the pro- 
posed methods meet these requirements, although considerable further 
advance can be made as better knowledge is obtained. For the condition 
of transporting water containing heavy sediment load, the reliability of 
the method is less certain, due to lack of knowledge of the reliability of 
formulae for the transportation of sediment, but for very heavy sedi- 
ment loads they are believed to be more reliable than methods which do 
not contain the concentration of sediment as a factor. As knowledge of 
the science of sediment transportation increases, the reliability of the 
method in this part of the field will be increased. 

Advantages of the author’s method of analysis are (1) that it applies 
to all cases of stable channels, not only those in self bourne alluvium, 
(2) that it introduces no new emperical coefficients, but deals only in 
such well known concepts as shear values and angle of repose, (3) that 
it shows what information is necessary to further perfect the method, 
(4) that the refinements of new knowledge can be introduced without 
modifying other aspects of the method and (5) that recent unpublished 
analysis has shown that the principles developed can be used to explain 
the shape of cross section which many canals take as the water flows 
through them. 

The author wishes to acknowledge the assistance of Prof. P. N. Lin 
in preparing this closing discussion. 
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DISCUSSION OF UNSTEADY FLOW IN OPEN-TYPE 
PIPE IRRIGATION SYSTEMS 
PROCEEDINGS-SEPARATE NO. 369 


E. H. TAYLOR,! A.M. ASCE, A. F. PILLSBURY,2 T. O. ELLIS,3 
and G. A. BEKEY.4—The very able discussion of the writer’s paper by 
Mr. Glover is gots appreciated, although the apparent attempt, here 
and elsewhere,» to glorify the effect of covers, and to minimize the 
effect of vents is a'bit puzzling. Actually, those parts of the Coachella 
system where surging was a problem were made operable by a com- 
bination of covers and of gates through the baffles. No mention was 
made of the operating difficulties with covers, to-wit: (a) surge some- 
times occurring at certain flows regardless of covers; (b) the blowing 
off of the covers from time to time; or (c) the necessity of decreasing the 
effectiveness of a given cover in order to obtain satisfactory delivery 
at that point. As for vents, model experiments would confirm that such 
field trials as were undertaken would not prove satisfactory. Until there 
are adequate field trials, the effectiveness of prototype vents must re- 
main pure speculation. 


1. Associate Prof. of Eng., Univ. of California, Los Angeles, Calif. 

. Prof. of Irrigation, Univ. of California, Los Angeles, Calif. 

. Engr. Assistant. 

. Junior Engr. 

. C. S. Hale, P. W. Terrell, R. E. Glover, and W. P. Simmons, Jr. 
Surge Control of the Coachella Pipe Distribution System. U.S. 
Bureau of Reclamation, Denver, Colo. Engineering monograph 
No. 17, Jan., 1954. 

. C. S. Hale, R. E. Glover, P. W. Terrell, and W. P. Simmons, Jr. 
Control of Surging in Concrete Pipe Distribution Systems. A.C.I. 
Proc. title No. 50-33, Mar., 1954. 
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DISCUSSION OF HUMID AREA SOILS AND MOISTURE FACTORS 
FOR IRRIGATION DESIGN 
PROCEEDINGS-SEPARATE NO. 426 


GEORGE H. HARGREAVEs!, A.M. ASCE--There has long been a need 
for a better understanding of soil, plant, and water relationships and 
their use in irrigation design. The engineer designing a canal system 
or a sprinkler system must know the seasonal amounts of water to be 
delivered and the required rates of delivery for specified periods. Mr. 
Larson has clearly outlined the crop, soil and moisture factors neces- 
sary for proper design. In the use of these factors, however, not only 
is caution required because of local deviations from average conditions 
but also because of the limitations of some of the factors themselves. 


Moisture Level and Yield—A large number of experiments show that for 
most crops growth is largely independent of the amount of moisture as 
long as there is readily available moisture in the crop root zone. Plants 
appear to take moisture as readily when the moisture content is only 
slightly above the permanent wilting point as when it is at the field ca- 
pacity. Veihmeyer and Hendrickson? studied the effect of moisture 
level upon the yield and quality of crops and found no significant differ- 
ence until the soil moisture in contact with the roots approached the 
permanent wilting point. Potatoes are known to be an outstanding ex- 
ception to the above rule as well as certain other coarse rooted plants 
such as celery, corn and lettuce. It seems probable that although irri- 
gation applications for some crops should be made when the soil mois- 
ture level is fairly high, for many crops this practice would lead to a 
waste of labor and water. 


Moisture Availability from Soils—Soil texture is an important factor 
determining the amount of moisture readily available to the crops. 

There also seems to be a definite relationship between the origin and 
physical and chemical properties of soils and their moisture relation- 
ships (ability to hold readily available moisture). Soils mapped as sand 
have been known to hold as much as 1.4 inches of readily available mois- 
ture per foot of depth. Medium textured soils (loams and clay loams) 
hold from 0.7 to 2.7 inches per foot while some clays hold as much as 
3.0 inches. Because of these wide variations average values based upon 
soil texture should be used with reservation. 


1. Civ. Engr., Inst. of Inter-American Affairs, Port-au-Prince, Haiti. 

2. Veihmeyer, F. J. and Hendrickson, A. H. “Water Holding Capacity 
of Soils, and its Effect on Irrigation Practices”, Agricultural Engi- 
neering, Vol. 19, No. 11, 1938 and “Irrigation Experiments with 
Grapes”, Calif. Agr. Exp. Station, Bull. 728, 1950. 
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TABLE XI 
MONTHLY AND SEASONAL CONSUMPTIVE-USE COEFFICIENTS 


Davis, Californie 


MONTHLY CONSUMPTIVE USE COEFFICIENTS "k" : SEASONAL 
Mer.: Apr.: May : June: July: Aug.t Sept.: Oct.: Now. : "K" 


CROP 


0.37: 0.753 0,81: 0.90: 1,04: 0,93: 0,64 
0.15: O,11: 0,15: 0,25: 0,65: 1.12: 


: 
Alfalfa 0.41: 0.33 0.76 


Asparagus 1.01: 0,41 0.57 


Fy 
Beans 0.21: 0,39: 0,91: 0,54 


n 


Cantaloupes 0.41: 0,52: 0,54: 0.52 
: 


Carrots 0,69: 0,91: 0,59: 0,49 
0.41 


0.37 


: 
Celery 0.21: 0,19: 0,35: 


Cora 


3 
0.17: 0,54: 0,58: 
3 


0.99: 1.00: 0,76: 0.67 


duous) 


Grain & Hay 


2 


0.52 


: 
0,17: 
: 


Grepes (Muscat) 0,352 0.31 


3 
Ladino Clover 1.03: 0,90 


Onions (early) 0,38 


Cnions (late) 0.39 


0,43 


: 
3 
: 

Fruit (Deci- : 
FY 
3 

Peas 

3 


Potstoes (early 0.77 
: 
Tomatoes 0,44: 0,98: 0.88 
: 


0,73 


a) Besed on consumptive use date published in "Suggested Subject Matter for Presenta- 
tion at Irrigotion Meetings" by L. J. Booher, 1948, Col. of Agr., Ext. Serv., 
Davis, California. File 12.1. 


: 7 
: : 
: 
: $ { 
4 
: 
: 
: 
: 
: 3 q 
: 0.34 : 0,10: || 
: 
0.47: 0.61: 
: : : { 
: : $ 3 3 
: 
: 
3 { 
3 4 
: : : : : : 
3 0.24: 0,38: 0,41: 0,41: 3 
: : : : : 
2 q 
3 i 5 
: 
| 
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Consumptive Use Coefficients—Experiment stations need to make many 
more measurements of consumptive use. Due to the present lack of 


measurements it is frequently necessary to apply data collected in one 
area to another. In order to do this successfully considerable knowledge 
of the factors which influence consumptive-use is necessary. Consump- 
tive-use for a specific crop varies with the rate and stage of growth of 
the plant and its leaf surface area as well as with the climatic factors 

of temperature, relative humidity, wind and daylight hours. Because 

all of these factors vary considerably throughout the growing season, 
irrigation design should be based upon monthly or short period values 

of consumptive-use rather than upon seasonal consumptive use coeffi- 
cients. 

In the computation of monthly rates of consumptive-use at a given 
location based upon known data from another area several methods may 
be employed. Where U. S. Weather Bureau pan evaporation is available 
at both locations a comparison of evaporation rates forms a good index 
of relative consumptive-use rates. Where evaporation data are not 
available either the Lowery-Johnson® or the Blaney-Criddle methods 
may be used. The Lowery-Johnson method offers some advantages 
where differences in relative humidity are of importance. Maximum 
daily temperatures (employed in the Lowery-Johnson method) and con- 
sumptive use both increase as relative humidity decreases. The Blaney- 
Criddle method offers the advantages of ease of computation and availa- 
bility of data. Unfortunately monthly consumptive-use coefficients are 
available for only a few crops and in only a few locations. 

The use of seasonal consumptive-use coefficients, given by Mr. 
Larson in Table IV, does not necessarily result in consumptive-use 
values that correlate closely with monthly or peak consumptive-use re- 
quirements to be used in irrigation design. Table XI gives monthly and 
seasonal consumptive-use coefficients computed for Davis, California 
by using the Bianey-Criddle method. Consumptive-use coefficients for 
months of maximum growth and leaf surface area average 48 per cent 
higher than seasonal coefficients. This difference, although of consid- 
erable importance, is reduced in humid areas. 

HARRY F. BLANEY*»°, M. ASCE—This paper provides a long needed 
stimulus for additional research in irrigation requirements in humid 
areas. This is a matter of importance not only in eastern United States 
but in areas throughout the world where irrigated agricultural devel- 
opment'is being expanded. Mr. Larson mentions the Blaney-Criddle 
method of determining consumptive use of water from climatological 
data in the text and in tables III, IV, and V, but leaves the reader in 
doubt as to the procedure. (1) 

Briefly, the procedure is to correlate existing consumptive-use data 
with monthly temperature, percentage of daytime hours, precipitation, 
frost-free (growing) period, or irrigation season. The coefficients so 


3. ASCE Transactions, Vol. 107, 1942, Page 1243. 
4. Principal Irrigation Engr., Western Section of Soil and Water Manage- 
ment, Soil and Water Conservation Research Branch, Agricultural 
Research Service, U. S. Dept. of Agriculture, Los Angeles, Calif. 

5. Member, Executive Committee, Irrigation and Drainage Div., ASCE. 
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developed for different crops are used to translocate or transpose con- 
sumptive-use data from one section to other areas in which climatologi- 
cal data alone are available. Research studies indicate that consump- 
tive use varies with the temperature and the extent of daytime hours, 
and with the available moisture (precipitation, irrigation, and ground 
water). By multiplying the mean monthly temperature (t) by the monthly 
percentage of daytime hours of the year (p), a monthly consumptive-use 
factor (f) is obtained. Then it is assumed that the consumptive use 
varies directly as this factor when an ample water supply is available. 
Expressed mathematically, U= KF = £ kf, in which U is the consump- 
tive use of crop in inches for any period; F is the sum of the monthly 
consumptive-use factors for the period (sum of the products of mean 
monthly temperature and monthly percentage of daytime hours of the 
year); K is the empirical consumptive-use coefficient (irrigation season 
or growing period); and the other terms are as previously defined. 

In table IV Mr. Larson shows tentative seasonal coefficients “K” for 
formula (U = KF = Consumptive Use in inches). In 1950 the writer in- 
spected the irrigated areas along the Atlantic Coast from Florida to 
Massachusetts, and suggested tentative values of “K” in connection with 
the design of sprinkler irrigation systems. These coefficients are in 
agreement in most instances with those shown in table IV. 

An analysis of the available data in the South Atlantic area indicates 
that the consumptive-use coefficients should be greater in the inland 
areas than along the coast. Table 1 shows the tentative coefficients 
proposed by the writer for the South Atlantic area. Additional research 
studies are needed to develop monthly cvefficients (k) for various crops 
grown in the Southern states. Table 2 illustrates the method of comput- 
ing the monthly consumptive use and irrigation requirements when data 
are available. 
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logical and Irrigation Data”, by Harry F. Blaney and Wayne D. Criddle, 
Div. of Irrig. and Water Conservation, Soil Conservation Service, 

U. S. Dept. of Agriculture, SCS-TP-96, Washington, D. C., 1950. 


q 
ing 
q 
| 
522-28 


Table 1. Tentative seasonal coefficients K for use in the consumptive- 
use formula U = KF, in the South Atlantic area of the United 
States. 


Period 


Winter 
Frost-free 
Corn 


Cot ton 


Deciduous (Orchard) 


Grapefruit 
Grain (small) Winter 
Frost-free 
Grass pasture Winter 
Summer 
Potatoes Frost—free 
Rice 
Strawberries Wint er 
Frost-free 
" " 


Tomatoes 


Vegetables 


Watermelons 8 


| 
Crop > : Coastal : Inland 
: (K) : _(K) 
0.40 0.60 
055 
-60 265 
60 | 
-80 
265 
1.00 
250 250 
055 
.60 
" " 255 
055 
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Table 2. Computed monthly consumptive use and irrigation requirement 
for grass pasture, Charleston, S. C., for dry year 1925. 


* Mean Monthly Monthly Monthly ‘Monthly “Consump-' Irriga- 


percent consumptive-' consumptive- precipi-. consump- tive use, tion 

Month? temper-' of use factor use * tation ztive use’ minus ;Fequire- 
hours: (£) (r) & (u-r) } (i) 


oF. Inches Inches 
59.2 8.36 J 1.28 2.97 1.69 264 
66.8 8.77 1.89 3.51 1.62 2.3 
71.0 9.67 1.96 
79.8 9.63 ° 5.49 4.99 

July 82.8 9.83 ° 2.38 5-70 

Aug. 81.2 9.31 ° 1.62 5.29 

Sept. 77.0 8.34 1.94 

Oct. 68.8 7.91 ‘ 3.08 3.26 


Season 
3/1 - 10/1 Total 34035 


7/ Based on field-irrigation efficiency of 70 percent. 
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